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A symmetric encryption method based on properties of quasicrystals is proposed. The ad- 
vantages of the cipher are strict aperiodicity and everywhere discontinuous property as well 
as the speed of computation, simplicity of implementation and a straightforward possibility 
of extending the method to encryption of higher dimensional data. 
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1 Introduction 

, Here we address the 'classical' cryptographic problem, i.e. protection of transmitted and stored 

data from divulgence and distortion. Such problem frequently arise when data creation or/and 
O . data usage are disjoint in time or/and space. There exist many approaches to this problem in 

modern cryptography, which can be divided into three main groups in accordance with their key 
primitives, i.e., unkeyed, symmetryc-key and public-key encodings, see [2ll31[9]. 

In the present work we propose a new generic encoding procedure based on the aperiodic 
CS) '■ point sets, called quasicrystals in the physics literature and model sets in the mathematics 

I literature [H |3] . Generally speaking, such an algorithm is a symmetric stream cipher endowed 

with strict aperiodicity (no periodic subsets) and everywhere discontinuous properties. 

The paper starts by preliminaries and general encoding idea, which is followed by the precise 
I statement of pertinent mathematical ingredients. Sections II and III. In Section [J] applications of 

O ' the proposed cipher to encoding of bitmap pictures are presented and discussed. The advantages 

of the approach are in the speed of computation and a straightforward possibility of extending 
^ \ the method to encryption of 3 and higher dimensional data. Possible range of private keys and 

^ ' some other applications of quasicrystals are discussed in Conclusions. 

2 Preliminaries and encryption procedure 

We call cut-and-project quasicrystal a discrete deterministic aperiodic point set A, in a finite- 
dimensional real Euclidean space. In many ways such quasicrystals resemble lattices in all 
but the translation symmetry. Sometimes they are even called aperiodic lattices or aperiodic 
crystals. In this work we make use of two remarkable properties of the quasicrystals: (i) no 
periodic subsets of any kind are contained in a quasicrystal, and (ii) the discontinuity of the 
'star map' between a quasicrystal A and its 'acceptance window' Q. 

General idea of the proposed approach is to take given digital data, set a one-to-one corre- 
spondence between information bites and integer numbers N, and then to map points of on 
a quasicrystal fragment A (see the first two lines of Fig. [1]). At this stage the data is mapped 
to the quasicrystal points. It is followed by the application of the star map to the quasicrystal 
points (see the second and the third lines of Fig. [1]), which is a crucial step of the encryption. As 
a result, we get the quasicrystal points, carrying our data, 'tossed up' in the acceptance window 
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J7. It remains to map the content of on a desirable lattice for the output (see the last two 
lines of Fig. [1]). 

Decryption proceeds in the opposite direction. 

The set of private keys in our encryption method is infinite and is discussed in Conclusions. 
In addition, one may choose to use any common substitution scheme. 

Generation of quasicrystal points is very flexible. It can either be done as a single stream or 
split into several blocks. Moreover, calculation of the points of each block is easily amenable to 
parallel computing. 

Using quasicrystals, all calculations can be carried out in integers. That guarantees absolute 
accuracy of coding and decoding procedures, assures the stability of the algorithm. 
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Figure 1: Scheme of the simplest quasicrystal-based cipher 



3 Mathematical ingredients 
3.1 Theoretical background 

Theoretical background of the work can be found in [U O O [71 [HI [TO]. 

Our aim is to consider two-dimensional data. We start with a choice of two straight lines 
Vi: y = Tx and V2 : y = t'x on a square lattice Z2 (see Fig. [2]), where r and r' are two irrational 
numbers to be specified later. The simplest choice is r = i(l + VE) and r' = i(l - Vb). One 
of the lines, say Vi, plays the role of the quasicrystal space onto which a 'cut' of the lattice Z2 
is projected. In the other direction V2, we choose a finite interval 0. Point of the lattice Z2 
becomes a point of the quasicrystal, after projection on Vi, provided its projection on V2 falls 
within fi, i.e. it is within the 'cut'. 

We denote by Z[t] all the real numbers of the form (a + br), with integers a and b. 

The star map between a quasicrystal point (a + br) and the point (a + br') of J7 is given by 

* : Z[r] — > Z[r'], 
(a + br) {a + br'), a,bGZ. 

There is an unlimited number of choices of the irrational pairs to use in our construction. 
For example, one may choose as r and r' the solutions of the two infinite series of the quadratic 
equations = mx + 1, m = 1, 2, . . . and = mx — 1, m = 3, 4, . . .. 

The one-dimensional quasicrystal or cut and project set At-(J7) is described as follows, 

Definition 1. Let be a finite interval and r, r' be irrational numbers, then 

Ar{Q) = {a + bT \ a,b e Z, a + br' e 
where is called the acceptance window of A,- ($7). 
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Figure 2: Construction of one-dimensional quasicrystal 

Let us underline, that At-{^) defined here is an infinite, uniformly dense, and uniformly 
discrete point set, while its window Q is a finite interval densely covered by the star map of the 
points of At- ($7). 

From the definition of AT-(fi) a number of properties of A,- ($7) can be shown [3l [10]. In 
particular, the aperiodicity, discontinuity of the star map and the existence of only two or three 
different distances between adjacent points of At-(Q). More precisely, we recall the following 
theorems. 

Theorem 1. Let 0, be a semi-closed interval. For every A,-(f7) there exist positive num- 
bers Ai,A2 € '^[t], such that the distances between adjacent points in At-(J7) take values in 
{Ai, A2, Ai + A2}. The distances depend only on t,t' and length of the acceptance window. 

More can be said about the distances between adjacent points, namely 



' {Ai, A2, A1 + A2}, 



{Xn+l - Xn\n eZ} = < 



when A 



A2 > 



{Ai, A2}, 
when A* — A 



One has a general scaling property of quasicrystals: 

Theorem 2. For every irrational numbers t,t', t ^ t' and a bounded interval ^l, there exist 
f' G (—1,0), f > 0, s G M and 17 (satisfying max{l + f', —f'} < |$7| <1), such that 

Ar{n) = sAf{n). 



Moreover, the distances between adjacent points in Af(il) are equal to f, 1 + f, 1 + 2f, when 
\Cl\ / 1 and to f, 1 + f, when = 1. 
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Finally we underline that the 1-dimensional quasicrystal At^^I) is an infinite discrete aperiodic 
set of points uniquely determined by the choice of r, r' and of acceptance window $7 = [c,c + d), 
where c and d are any real numbers. Due to the theorem 2, without loss of generality, we can 
consider r' G (-1, 0), r > 0, c = and d e (max{l + r', -r'}, 1]. 

In higher dimensions, it is adequate for our purposes to take for the quasicrystals a straight- 
forward concatenation of one-dimensional ones in pairwise orthogonal directions. 

3.2 Computational aspects 

First let us introduce two specific examples of one-dimensional quasicrystals A,- (17). 

Example 1. Take r = , t' = , the two roots of the quadratic equation = x -|- 1. 
Then t + t' = 1, rr' = — 1. Choosing Q = [0,d), d = ^ the conditions of Theorem 2 are satisfied: 
r' « -0.62 E (-1,0) and r ^ 1.62 > 0. 

Due to Theorem 2, the distances between adjacent points of the quasicrystal AT-{[0,d)) are 
Ai = r, A2 = 1 T and A3 = 1 2r. 

As the seed point of the quasicrystal Ar{[0,d)), we can choose any point x = a + br, such 
that x' = a + br' G [0, d). For our example we put x = 0. 

Unlike Definition 1, which is not constructive, here we have the information for fast generation 
of quasicrystal points. We can move right or left from the seed point by adding or subtracting 
one- by-one of the distances Ai, A2 or A3. 

Suppose we have already established, that x is a point of At-([0, d)), i.e. x' G [0, d). The point 
adjacent to x from the right is one of the three 

X + T, X -|- 1 -|- r or X -|- 1 -|- 2r. 

In order to decide which one is the case, we verify one-by-one the inclusions x -|- r' E [0, d), 
X -|- 1 -|- r' G [0,d) or X -|- 1 -|- 2r' G [0, d). The first confirmed inclusion determines the new 
quasicrystal point. 

In such a way the following two subsets of the quasicrystal A^([0, d)) were obtained 



Here the second line of each example contains the corresponding sequence of distances between 
quasicrystal points. The two quasicrystals differ by the sequence in which the three distances 
Ai, A2 and A3 were tried in the construction. 

The first case is was shown as the second line on Fig. [TJ 

Example 2. Next quasicrystal AT-(r2) is built using the irrationalities t = 1 + \/2, r' = 1 — \/2, 
roots of the quadratic equation x^ = 2x -|- 1. 

The conditions of Theorem 2 are again satisfied: 

r' « -0.42 G (-1,0) and r 2.42 > 0. Putting 0, = [0,1), there are only two distances 
Ai = T and A2 = 1 -|- r between adjacent points of the quasicrystal At-([0, 1)). 

Using the zero point as the seed point, we obtain the cut and project set 



{ 



-l-2r, -T, 0, 1+T, 2+2t, 2+3t, 3+4t,...}^ 
^ {..., A2, Ai, A2, A2, Ai, A2, ...}; 

r, 0, 1 + r, 2 + 3r, 3 + 4t,...}^ 

^ {...,Ai, A2, A3, A2,...}. 



{• 



{. . . ,-l-3r,-2T,-T, 0, l+T, 1+2t, 2+3t, 2+4t, 3+5t, . . .} 
^{..., A2, Ai, Ai, A2, Ai, A2, Ai, A2, ...}. 
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It may happen that some point x' is arbitrary close to an end point of Jl. Then it would be 
a time demanding computational task to decide whether x' G or x' ^ fi. Such a difficulty is 
simply avoided by disqualifying from our consideration any point x' which comes closer to the 
boundary of than the distance e agreed in advance. 

The real numbers of the form a + br can be understood as given by two integer components 
(a, b). Since our operations are only addition or subtraction, the arithmetics of such numbers is 
elementary introduced. Consequently all transformations are performed with absolute precision. 

The transformations between the original and encrypted data go through the several stages. 
At each stage we need only the previous point. Therefore we have only minimal requirements 
to RAM. 

4 Application to image coding 

In this section fragments of quasicrystal point sets are used for encryption/decryption of 2- 
dimensional digital data. 

Assuming that the data is given on the rectangular grid containing final number ni, n2 
of points in each direction, we deal with final fragments A(*)(Qj), i = 1,2 of the quasicrystal 
A(f2) = A''"'^'*(r2i) (8) A*^^^(r22)- Number of points in each fragment is given by digital data size 
ni = |A»(Oi)l- 

Each is an interval on real axis in mutually orthogonal directions. They can be chosen for 
both fragments independently. 

Our main tool is the one-to-one map (star map) between the finite point sets Aj(r2j) and its 
image in Vti 

Ai(Oi) i — > Qi, i = l,2. 

Note, that roles of Aj(rij) and its image in fij can be interchanged, because the fragments 
Aj(Qj) are finite sets of points. 

Fast generation of Aj(J7j) was described in previous section. It requires that one is given the 
values of ni, n2, Tj, t| and fij, and two seed points of Aj(J7j), i = 1,2. 

Suppose one is given a bitmap picture P{L), sampled on the points {x, y) of a fragment L of 
a rectangular lattice. 

Then we treat independently each coordinate of the data using one fragment of the quasicrys- 
tal as illustrated in Fig. [H 

4.1 Basic encryption algorithm 

• Construction of two different 1-dimensional quasicrystals A^^^ and A^^^ of sizes ni and n2 
appropriate for the given data P{L). (Each fragment, as well as its acceptance window- 
interval can be situated anywhere on the corresponding real axes.) 

• Matching one-by-one the sequence of coordinates x with points of A'-^^ similarly matching 
the coordinates y with points of A^^^. (In this way x becomes point xi + X2T of the A^^^ 
and y becomes point yi + y2T of the A^^^ where xi, X2, yi, 2/2 are integers.) 

• The main step of the encoding is the star map of A^^^ A^^^ into the corresponding window 

fil fi2 

AW(Oi)^r?i, A(2)(j^2) -^5^2; 

xi + X2T1 ^ xi + X2t[, yi + y2T2 ^ yi + y2T2- 
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(The sequence of the original data coordinates is dramatically changed due to the discon- 
tinuity of the star map.) 

• It remains to map the points of encrypted quasicrystal fragments, which are now in 0,i 
and to a sequence of integers along x and y directions. 

The present algorithm permutes the data points without changing the value of the data at 
each point. It is straightforward to consider the possible values of the data function as another 
discrete point set and to apply the same algorithm to it. Effectively we are thus dealing with 
a three-dimensional problem, where the data function takes only two values at 3-dimensional 
points. 

All the maps we have used in the encryption algorithm are one-to-one, therefore they can 
be used in the inverse order for the decryption. For the same reason, nothing prevent us from 
iterative using our algorithm. 

In certain types of data the encryption algorithm leaves residue of the orthogonal directions, 
see Fig. [8)3. It can be avoided in a many different ways. One of them was mentioned above, when 
the problem was viewed as a three-dimensional one. Another possibility is illustrated on Fig. [3h. 
The two-dimensional data is taken as an ordered one-dimensional sequence (modification 2). 
The particularly fast possibility is to introduce the additional cyclic permutations of the points 
(modification 1), see Fig. [3]i. 




abed 



Figure 3: Application of the algorithms proposed in Section [Ha - initial image; b - main 
algorithm; c - second modification; d - first modification 

The computation time of our algorithms is compared in Table 1. Computations were carried 
out by Pentium M processor 1.70 GHz with 1 Gb RAM. As the test image we took the 512x512- 
pixel 8 bit grayscale picture Lenna (see Figs. |4]and[5l). But the second modification was tested 
on 200x200-pixel 8 bit grayscale picture Lenna (Fig. 5). Note, that the program does not cut 
the images into blocks for processing. 



Table 1. Program working time. 





Main algorithm 


Modification 1 


Modification 2 


time 


3 seconds 


6 seconds 


192 seconds 



Program realizing these algorithms can be freely downloaded from http://132.204-90.210/Tnaryna/ 
and listing of the main coding algorithm is adduced in Appendix lAl 

5 Conclusions 

In the paper one of the simplest versions of the quasicrystal-based symmetric cipher is described. 
Even in this form it appears adequate for many applications. Its complicated key is described 
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below. It is based on strongly aperiodic sets, includes everywhere discontinuous 'star map', and 
it is absolutely stable, since it operates with integer numbers only. Similarly, damaging a few 
data points does not affect the rest of the encryption/decryption process. 

Furthermore coordinates of the data can be divided into several non-overlapping intervals, 
then contents of all intervals could be encrypted in parallel. 

5.1 Cryptovariables 

Application of the algorithm requires that the following variables (cryptovariables of the private 
key) are fixed 

T jr' _aJ)_c-d-sml-X , where 

T and t' are irrational numbers, solutions of one of equations of the two infinite series of the 
quadratic equations x'^ = mx + 1, m = 1, 2, . . . or = mx — 1, m = 3, 4, . . .; 

a and b are integers, fixing a seed point of the quasicrystal; 

c and d are arbitrary real numbers, fixing the position of the acceptance window fi; 

X takes values {'+', '-'}, defining direction of quasicrystal construction relative to the seed 
point; 

sml sequence of distances of adjacent points during quasicrystal generation; 

£ is the maximal distance from the end points of Q at which the star map points are still 
considered; 

/ is the positive integer number specifying number of iterations of the certain quasicrystal cipher. 

More elaborate encryption system can be built by exploiting properties of 1 or 2-dimensional 
quasicrystals. Using the present algorithm in several iterations, one may choose new quasicrystal 
for each iteration, or new position of the quasicrystal fragment, or different size, position of the 
acceptance window, etc. 

Rich scaling properties of the quasicrystal could be involved into encryption scheme [3]. 

Genuinely two-dimensional quasicrystals can be used for the encryption, some examples of 
the possibilities are found in [12]. The number of distinct quasicrystals of this kind is unlimited. 

One can easily build one-dimensional quasicrystal by projecting cuts of n-dimensional latices 
on the direction of the quasicrystal. 

Intriguing appears to be the possibility to read the initial lattice data as an ordered stream of 
points and map them into single one-dimensional quasicrystal, perform the star map and reorder 
transform points into encrypted data. There are also other applications of the quasicrystals e.g. 
for generation of random numbers [6j, even some cryptographic ones pT] . 
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A Program listing 

Realization of the main algorithm on Delphi 6.0: 
unit Unitl; 
interface 

uses 

Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms, 
Dialogs, StdCtrls; 
type 

TNode = Record 
a : integer; 
b : integer; 
rs, rb : real; 
ns, nb : integer; 
end; 

TSequence = Record 

n : integer; 

Nodes : array of TNode; 

a,b : integer; 

s : string; 
end; 

TForml = class (TForm) 
Buttonl : TButton; 

Button2: TButton; 

procedure ButtonlClick(Sender: TObject); 
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procedure FormCreate(Sender: TObject); 
procedure Button2Click(Sender: TObject); 
private 

Private declarations 
public 

Public declarations 
end; 
var 

Forml : TForml; 

Roots, RootB : real; 
implementation 
{$R *.dfm} 

function CalcSmall (a, b: real) : real; 
begin 

CalcSmall := a + b * RootS; 
end; 

function CalcBig (a, b: real) : real; 
begin 

a := a + 1; 

b := - b; 

CalcBig := CalcSmall (a, b); 
CalcBig := a + b * RootB; 
end; 

function CheckShort (a, b: integer; var Node : Tnode; dir: integer) : boolean; 
var tempB : real; 

res : boolean; 
begin 

a := a - dir; 
b := b + dir; 
tempB := CalcBig (a,b); 
res := ((tempB < 0) and (tempB < 2)); 
If res then begin 
node. a := a; 
node.b := b; 
node.rb := tempB; 
end; 

CheckShort := res; 

end; //CheckShort (a, b: integer; var node : Tnode; dir: integer) : boolean; 
function ChcckMedium (a, b: integer; var node : Tnode; dir: integer) : boolean; 
var tempB : real; 

res : boolean; 
begin 

a := a + dir; 
tempB := CalcBig (a,b); 
res := ((tempB > 0) and (tempB < 2)); 
If res then begin 
node. a := a; 
node.b := b; 
node.rb := tempB; 
end; 
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CheckMedium := res; 
end; //CheckMedium (a, b: integer; var node : Tnodc; dir: integer) : boolean; 
function CheckLong (a, b: integer; var node : Tnode; dir: integer) : boolean; 
var tempB : real; 

res : boolean; 
begin 

b := b + dir; 
tempB := CalcBig (a,b); 
res := ((tempB > 0) and (tempB < 2)); 
If res then begin 
node. a := a; 
node.b := b; 
node.rb := tempB; 
end; 

CheckLong := res; 

end; //CheckLong (a, b: integer; var node : Tnode; dir: integer) : boolean; 
function FillSequence (var seq : TSequence; a, b, n, dir: integer): string; 
var res : string; 

aPrev, bPrev : integer; 

i : integer; 

begin 

res := "; 

seq.n := n; 

seq.a := a; 

seq.b := b; 

SetLength (seq.Nodes, n); 

seq.Nodes[0] 

seq.Nodes[0].b := b; 

seq.Nodes[0].rb := CalcBig (a, b); 

seq.Nodes[0].rs := CalcSmall (a, b); 

seq.Nodes[0].ns := 0; 

for i := 1 to n - 1 do begin 

aPrev := seq.Nodes [i-l].a; 

bPrev := seq.Nodes [i-l].b; 

seq.Nodes [i].ns := i; 

if dir = 1 then begin 

if CheckShort (aPrev, bPrev, seq.Nodes [i], dir) then 
res := res + 's' 

else if CheckMedium (aPrev, bPrev, seq.Nodes [i], dir) then 

res := res + 'm' 
else if CheckLong (aPrev, bPrev, seq.Nodes [i], dir) then 

res := res + '1' 

end else if CheckLong (aPrev, bPrev, seq.Nodes [i], dir) then 
res := res + '1' 

else if CheckMedium (aPrev, bPrev, seq.Nodes [i], dir) then 

res := res + 'm' 
else if CheckShort (aPrev, bPrev, seq.Nodes [i], dir) then 
res := res + 's'; 
seq.Nodes [i].rs := CalcSmall (seq.Nodes[i].a, seq.Nodes [i].b); 
end; 
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seq.s := res; 

FillSequence := res; 
end; //FillSequence (var seq : TSequence; a, b, n, dir: integer): string; 
Procedure SaveSequence (seq : TSequence; fname : string); 
var f : TextFile; 
var i : integer; 
begin 

AssignFile (f, fname); 

rewrite (f); 

writeln (f, seq.n,' : ', seq.s); 
for i:= to seq.n - 1 do 

writeln (f, seq.Nodes[i].ns, ', ', seq.Nodes[i].nb, ', ', seq.Nodes[i].a, ', ', seq.Nodes[i].b, ', 
', seq.Nodes[i].rs, ', ', seq.Nodes[i].rb); 
Flush (f); 
CloseFile (f); 

end; //SaveSequence (seq : TSequence; fname : string); 
Procedure SortSequence (var seq : TSequence; dir: integer); 
var node : TNode; 

i,j : integer; 

flag : integer; 
begin 

for i := seq.n-2 downto do 
for j:= to i do begin 
flag := 0; 

if dir = 1 then begin 

if seq.Nodes[j].rb > seq.Nodes[j + l].rb then 
flag := 1; 
end else 

if seq.NodespJ.rs > seq.Nodes[j + l].rs then 
flag := 1; 
if flag = 1 then begin 

node := seq.Nodes[j + 1]; 
seq.Nodesp + 1] := seq.Nodesp]; 
seq.Nodesp] := node; 
end; 
end; 
if dir = 1 then 

for i := to seq.n do 
seq.Nodes[i].nb := i; 
end; //SorteSequence (var seq : TSequence; dir: integer) 
Procedure Codelmg (dir: integer); 
var fnamel, fname2 : string; 
Bitmapl, Bitmap2 : TBitmap; 
SeqX, SeqY : TSequence; 

i,j : integer; 

CurWidth, CurHeight : integer; 
strX, strY : string; 

begin 

if dir = 1 then begin 
fnamel := 'in.bmp'; 
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fname2 := 'out.bmp'; 
end else begin 

fnamel := 'out.bmp'; 

fname2 := 'back.bmp'; 
end; 

Bitmapl := TBitmap. Create; 
Bitmapl.LoadFromFile (fnamel); 
Bitmap2 := TBitmap. Create; 
CurWidth := Bitmapl.Width; 
CurHcight := Bitmapl. Height; 
Bitmap2.Width := CurWidth; 
Bitmap2. Height := CurHeight; 
StrX := FiUSequence (SeqX, 1, 0, CurWidth, 1); 
StrY := FiUSequence (SeqY, 1, 0, CurHeight, 1); 
SortSequence (SeqX, 1); 
SortSequence (SeqY, 1); 
for i:= to CurWidth - 1 do 
for j:= to CurHeight - 1 do 
if dir = 1 then 

Bitmap2.Canvas.Pixels [i,j] := Bitmapl. Canvas.Pixels [SeqX.Nodes [i].ns, 

SeqY.Nodes p].ns] 

else 

Bitmap2. Canvas. Pixels [SeqX.Nodes [i].ns, SeqY.Nodes |j].ns] := Bitmapl. Canvas. Pixels 

[ij]; 

Bitmap2.SaveToFile (fname2); 

Bitmapl. Free; 

Bitmap2.Free; 

if dir = 1 then begin 

SortSequence (SeqX, -1); 

SortSequence (SeqY, -1); 

SaveSequence (SeqX, 'x.txt'); 

SaveSequence (SeqY, 'y.txt'); 
end; 

end; //Codelmg (dir: integer) 

procedure TForml.ButtonlClick(Sender: TObject); 
begin 

Buttonl. Enabled := False; 
Button2. Enabled := False; 
Codelmg (1); 
Buttonl. Enabled := True; 
Button2.Enabled := True; 
end; 

procedure TForml.FormCreate(Sender: TObject); 
begin 

Roots := (1 + sqrt (5)) / 2; 
RootB := (1 - sqrt (5)) / 2; 
end; 

procedure TForml.Button2Click(Sender: TObject); 
begin 

Buttonl. Enabled := False; 
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Button2. Enabled := False; 

Codelmg (-1); 

Buttonl. Enabled := True; 

Button2.Enabled := True; 
end; 
end. 

B Application examples 

In this section we will show several examples of the application of the main algorithm, its 

modification 1, and modification 2, as well as its iterations several times. Let us emphasize that 
the data was encrypted as a single block. Computing time is ranging between two seconds (one 
iteration) and nine seconds (five iterations) for the main algorithm. 




Figure 4: Main algorithm coding without iterations 




Figure 5: Modification 1 of the algorithm without iterations 
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Figure 6: Modification 2 of the algoritlim without iterations 



An omnibus across the bridge 

Crawls like a j^iiow butterfly, 

And.tiene and there, a passer-by 
Shows like a little restless midge. 



Big tSarges full of yelbw hoy 

Aft mooKd 

against the shadow 

And like a TysXXoss silken scarf, 
The thick fog hangs along the quay. 

The yellow leaves begin to fade 
And flutter from the Temple elms 
And at my feet the pale greenThames 
Lies like a rod of rippled jade. 




- 

'l.:.l£ 



... 


. ■■■■■■ . .■;--v.-a^ 





















Figure 7: Main algorithm coding without iterations 



An omnibus across the bridge 

Crawls like a yellow butterfly, 

And, here and there, a passer-by 
Shows like a little restless midge. 



Big barges full of yellow hoy 

Are moored 

against the shadow 

wharf 

And like a \\ieXXom silken scarf, 
The thick fog hangs along the quay. 

The yellow leaves begin to fade 

And flutter from the Temple elms 
And at my feet the pale green Thames 
Lies like a rod of rippled jade. 

Figure 8: Modification 1 of the algorithm without iterations 
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